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We consider the influence of mode - mode coupling in the inflaton field on the spectrum of 
primordial fluctuations. To this end, we formulate a phenomenological model where the inflaton 
fluctuations are treated as a fluid undergoing turbulent motion. Under suitable assumptions, it 
is possible to estimate the size and scale of fluctuations in velocity, which upon reheating induce 
corresponding fluctuations in the radiation energy density. For De Sitter inflation the resulting 
spectrum is scale invariant on all scales of interest. The amplitude of the resulting spectrum is 
compatible with known observational limits. This suggests that the hypothesis of a extremely 
weakly coupled inflation could be relaxed without affecting the predictions of the model. Principal 
PACS No 98.80.Cq; additional PACS nos: 04.62.+v, ll.10.Wx, 47.27 
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The object of this paper is to present models of fully nonlinear fluctuations during the inflationary era conducing 
. to a scale invariant primordial density contrast spectrum with an amplitude consistent with COBE observations. 
' Inflationary models were originally introduced as a solution for the so-called puzzles of standard hot Big Bang 
cosmology [Q, namely the horizon, flatness and photon to baryon ratio problem A period of exponential 
expansion of the Universe was postulated, whereby the ratio of the Universe increased by a factor e 60 or so, followed 
by a period of reheating, during which the temperature of radiation was raised up to a final value of 10 9 GeU or more 
(we assume units with K = c = k = 1) Q. From then on, cosmic evolution followed the lines of standard cosmology 

. Although many implementations of inflation have been proposed, most attention has been devoted to the simplest 
scenario, which was Linde's chaotic inflation with a single inflaton field ||. In this model, inflation is powered by a 
t^- . scalar field <f>, the inflaton, slowly rolling down an effective potential V (</>) . We shall simplify the model further by 
ON ■ assuming the Universe is described by a spatially flat, Friedmann - Robertson - Walker geometry. Inflation begins 
when the effective potential becomes the dominant form of energy density, and ends with the decay of the inflaton 
O onto radiation. During inflation the potential energy V (</>) acts as an effective cosmological constant. It can be said 
that no satisfactory contending explanation of the cosmic puzzles is available H . 

Soon after the original proposal, it was realized that inflation could perform a subtler task: to provide a framework 
for explaining the origin of primordial density fluctuations p|. Quantum fluctuations of the inflaton field distort the 
reheating surface, inducing a primordial density contrast (see Q, we shall review this argument in greater detail 
below) 

><: fc~§* a) 



H being the inflationary Hubble parameter. A similar argument shows that gravitational waves are also being created, 
with an amplitude h ~ Hlm v . The validity of Eq.(Q) has been corroborated by several different methods Jj| (for a 
dissenting viewpoint, see [Ell). 

In order to obtain a concrete prediction from Eq.(|l]) we must estimate the quantum fluctuations 5(f). The usual 
approach treats these fluctuations as a free field (for example, in the seminal paper by Starobinsky pd|]) in its (De 
Sitter invariant) vacuum state or very close to it JO . Then a simple calculation within quantum field theory in 
curved spaces allows us to evaluate the Scj) in Eq.(P~as 5<f) ~ H p3|. All quantities in the right hand side of Eq.(Q) 
are evaluated as the relevant mode leaves the horizon. There are well defined ways to relate the spectrum of scalar 
and tensor primordial fluctuations to a given potential, which are generally described as the potential reconstruction 
program | ju| . 

As our understanding of the formalism and the complexity of observational data progresses fl5f| , it becomes clear 
that simple potentials such as single powers of </> are not rich enough to account for observations. Thus the potential 
reconstruction program generally assumes more complex functional forms, depending upon several parameters fi~4|| . In 
these more general potentials their higher derivatives may not be negligible, leading to nonlinear interactions between 
fluctuations. In such a case, the usual way of estimating the primordial density contrast would be invalid. The 
relevance of non linear fluctuations was already taken up in Ref. Also several authors Jl7|-p0| have considered 

models where fluctuation self interactions were treated perturbatively. 
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The main goal of this paper is to develop ways of estimating the primordial density contrast in chaotic inflation 
models without presupposing that couplings among fluctuations are negligible. Of course, working out the structure 
of quantum fluctuations on even simple states of fully nonlinear quantum field theory is a daunting task |pl[ . However, 
a simpler alternative is available, namely, the application of hydrodynamics to describe the macroscopic behavior of 
quantum fluctuations. This is possible because these fluctuations, as far as it is relevant to our present concerns, may 
be described by a c-number energy momentum tensor subject both to the usual conservation laws and the Second 
Law of thermodynamics Jl3| . There is therefore an equivalent fluid description, consisting of a classical fluid whose 
energy momentum tensor and equation of state reproduce the observed ones for the quantum fluctuations. Solving 
the dynamics of this equivalent fluid yields answers to all relevant questions concerning the behavior of the actual 
quantum fluctuations. In this paper we shall not press the issue of whether the equivalent fluid is anything more than 
a convenient computational device. 

Along with the assumption of free inflaton fluctuations, we must question whether these fluctuations arc in their 
vacuum state. This assumption usually rests on the so called cosmic no hair theorem, which states that any reasonable 
initial quantum state for the field will quickly relax to a De Sitter invariant state, or its closest available analogue 
p2| . The application of this theorem in models with a hundred or more e-foldings of inflation is justified, but these 
models must be rejected on the grounds that they also predict a value of the density parameter f2 unacceptably 
close to one (see Appendix and Ref. [^3); it is also possible to develop inflationary models in open universes, see Ref. 
pi|). When the duration of inflation is close to minimal, it becomes likely that some large but observable scales will 
leave the horizon before the cosmic no hair theorem is able to operate. It is important therefore to set a realistic 
initial condition such as a Planckian distribution of particles with temperature T. The relevance of thermal effects in 
inflation was already pointed out (in a different context) in ref. [ p5| . 

An immediate consequence of energy momentum conservation and the Second Law is that when velocities are low, 
the phenomenological fluid may be described within the Eckart theory of dissipative fluids [^6[^7) (for an analysis of 
the limitations of Eckart 's theory see ^8|). It follows that it obeys a continuity equation and a curved space time 
Navier - Stokes equation. The model is then defined by giving the equation of state and the viscosity of the equivalent 
fluid. The advantage we gain is that these are features that can be computed locally. As far as the relevant scales 
are much below the curvature radius, it is possible to use for them their standard flat space time values. At high 
temperatures we obtain the equation of state for radiation, p = (l/3)p, and a dynamic viscosity r/ ~ T 3 . Since the 
speed of sound, being close to light speed, is much higher than the characteristic speed of the fluid, the flow may be 
considered incompressible. There will be fluctuations in velocity, nevertheless, and these are the ones responsible for 
density fluctuations, as it is well known 

As we shall show below, conditions in the early stages of inflation are such that, for generic initial conditions, the flow 
of the equivalent fluid is highly turbulent, meaning that the corresponding Reynolds number is well over a thousand. 
The possible role of turbulence in the formation of cosmic structures has been studied in some detail even before 
inflationary models were introduced (see |2£||3(J for a textbook account). These early attempts were abandoned 
because there were no natural mechanisms for the production of primordial turbulence, and the density contrast 
predicted was generally too high. After inflation was proposed, the matter was taken up again by Goldman and Canuto 
who studied longitudinal turbulence excited by density fluctuations in the radiation and matter dominated eras. 
Our work should not be understood as a continuation of this line of research, but rather as a variation on conventional 
models of primordial fluctuation generation (cfr. Jl|,^,|llJ) whereby quantum fluctuations of the inflaton field leave their 
imprint on the primordial density field and are subsequently washed away. In other words, the underlying physics in 
our model is the same as in these more familiar approaches to primordial fluctuation generation: we do not question 
the ultimate quantum origin of the fluctuations (another difference with the cosmic turbulence theory from the fifties), 
but simply borrow insights from hydrodynamics to describe the macroscopic behavior of these fluctuations, rather 
than rely on possibly oversimplified linearized microscopic models. The conditions of validity of our procedure are the 
assumptions that the energy momentum tensor of fluctuations is a c-number quantity (which ought to be true at any 
scale below Planck's) and the Second Law (which, unlike the Third, has not been challenged yet to our knowledge). 

A second goal of this paper is to demonstrate the application of the hydrodynamic approach by discussing some 
simple solutions to the non linear Navier - Stokes equation, and the resulting spectra of primordial fluctuations. Since 
the general solution to the Navier- Stokes equation is certainly not available, this requires the appeal to physical 
insights to simplify the problem. 

The basic mechanism of non linear hydrodynamic evolution is the energy transport between eddies of different size 
through mode - mode coupling, and the viscous dissipation of small scale eddies. The fundamental issue regarding 
model building is to obtain a closed form energy balance equation, which allows us to follow in time the shape of 
the energy spectrum. This usually involves some closure hypothesis to reduce the infinite hierarchy of dynamical 
equations for velocity correlation functions to a manageable set Lacking anything better, we shall fall back on 
the time-honored hypothesis that the effect of smaller eddies at any given scale may be simulated by a scale dependent 
effective viscosity 33-3f|; for concreteness, we shall follow Heisenberg's 1948 formulation of this idea |36j| . 
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The Heisenberg theory admits some very simple solutions with the property of self-similarity. These have been 
worked out by Chandrasekhar |5?J and generalized to Friedmann - Robertson - Walker (FRW) backgrounds by 
Tomita et al |Q. These solutions agree with the Kolmogorov 1941 theory in the inertial range failing to 
reproduce observations for very small eddies. Fortunately, we shall only require the solution in the opposite limit of 
very large eddies, where it is trustworthy^. With minor adjustments, Tomita's analysis of turbulence decay in FRW 
space times also provides a solution to the evolution of our equivalent fluid. Self-similarity is an appealing feature to 
us, as we do expect that a generic flow will be eventually brought to some sort of steady state by the inflationary 
expansion. 

The task at hand is then to study the evolution of a typical eddy as it is blown up by the universal expansion, 
exchanging and dissipating energy while inside the horizon, and freezing when outside, until it reaches the reheating 
hypersurface and delivers its energy to radiation. By assuming that the turbulent velocity fluctuations in the eddy 
produce fluctuations in the energy density of radiation in the usual way, we shall be able to relate the primordial 
density contrast to the features of the original self similar turbulence. The resulting spectrum may be matched 
against the known data on the cosmic microwave background ]l5|] , providing a crucial test of the inner consistency 
and viability of the approach. Our conclusion shall be that, insofar as the horizon remains constant during inflation, 
the spectrum of primordial density fluctuations produced by self similar flows is strictly scale invariant (n = 1, see 
friof ) at large scales. Quantitative agreement with observations may be obtained without any special fine tuning. 

A few comments on the possibility of deriving present day spectra from early Universe hydrodynamics are in order; 
after all, one of the main criticisms against the conventional cosmic turbulence theory was that, even if turbulence were 
efficiently generated in the Early Universe, it would decay and become uninteresting well before recombination. This 
criticism does not apply to the present work, since we shall show that the results of customary inflationary models and 
self similar flows concerning the primordial density field immediately after reheating are identical. Further evolution 
of the primordial density contrast after reheating is well described by the theory of linear density fluctuations in 
an expanding Universe, a subject properly covered by many textbooks (cfr. fl]|i| J^ , f4l[ ) ; the result, for both the 
conventional inflationary models and the present ones, is that these fluctuations in the primordial density contrast 
survive to recombination time. 

The key point in our analysis is the evolution of the flow during inflation, as the different modes leave the horizon 
and " freeze" , and we discuss this issue at some length. The amount of turbulence at reheating is ample enough to 
seed fluctuations at the level suggested by COBE data, provided the initial temperature is large enough. The lower 
bound in the temperature, however, is not so large that would invalidate the vacuum dominance condition which is 
a presupposition of Inflation. Of course, whether Inflation is likely to happen or not is a difficult question (cfr. [fl2[) 
lying beyond the scope of this paper. 

To conclude, the main objective of this paper is to show that it is possible to develop sensible models of inflation 
where inflaton fluctuations evolve nonlinearly and are very far from their vacuum states. The connection of the physics 
of primordial fluctuations to hydrodynamics opens up a wealth of new interesting phenomena, such as intermittency 
in the primordial spectrum |34|,^| and Burgers turbulence [ pf4| ] , with a strong potential impact on our understanding 
of the evolution of cosmic structures. Moreover, it is appealing to be able to account for a macroscopic phenomenon, 
such as fluctuation generation on super horizon scales, mostly on macroscopic terms (for an independent attempt 
in this direction, see [^5| ). Most importantly, by not assuming that higher derivatives of the potential are a priori 
negligible, we avoid a possible conflict within the potential reconstruction program |OL 

The rest of the paper is organized as follows. In next section we provide a brief summary of hydrodynamics in 
flat and expanding universes, in order to set up the language for the rest of the paper. In section III we proceed to 
discuss the equivalent fluid description of inflaton fluctuations, and how to extract the primordial density contrast 
therefrom. As a simple application of the method, we consider briefly the case of free fluctuations, showing that the 
model leads back to the conventional results. In section IV we present Chandrasekhar's self-similar solutions and their 
generalization to expanding universes, deriving the corresponding scale invariant primordial contrast; we also discuss 
how relevant these solutions are as actual descriptions of the inflationary period. We state our main conclusions in 
the final section. 



II. HYDRODYNAMIC FLOWS 



x We wish to point out that the applicability of Kolmogorov's spectrum to large scale turbulence should not be taken for 
granted Q 
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A. Flows in flat space time 



The equations governing the dynamics of a fluid in local thermodynamic equilibrium are the continuity and Navier- 
Stokes ones, which, in the case of flat space time, read: 

^ + (U.V)p = (2) 



^ + (U-V)U = --Vp + vV 2 XJ (3) 
ot p 

where we have assumed incompressibility, valid when typical velocities are much smaller than the sound velocity; 
v = r\j p is the kinematic shear viscosity. The transition from laminar to turbulent motion can be universally described 
by the dimensionless "Reynolds" number: 

R=^ 
v 

where U is a typical velocity and L a typical length scale. This number represents the order of magnitude of the ratio 
of the incrtial to the viscous term. Low Reynolds numbers correspond to laminar motion, while high ones suggest 
turbulent behavior. 

In general, the velocity profile displays variations in space and time. This implies that the flow must be described 
probabilistically. Thus, each quantity involved in is divided in its mean value and a fluctuation from it; for 

example, we write U = U + u, where u stands for the fluctuating part of the velocity. In the case where motion is 
isotropic, the mean value U for the velocity must be zero, since otherwise there would be a preferred direction. 

To analyze the system's behavior, we define the two-point one-time correlation function for the velocity: 

R ij (x,x',t) = (mixfiujix't)) (4) 

In the case of homogeneous and isotropic motion, this correlation must be only a function of the time t and the 
distance between x and x', i.e. Rij(x,x',t) = Rij(r,t), where r = |x — x'| . Observe that Ru(0,t) (summation over 
repeated indices must be understood) is twice the average energy density of the flow at time t. From (^) we obtain 
the equation that this correlation must obey, namely: 



0_ 

at 

where 



R tJ (r, t) = T tj (r, t) + P lJ (r, t) + 2uV 2 R tJ (r, t) (5) 



Pij(r,t) = - ( A (p(x,t) Uj (x',t)) - ^- (p(x',t) Ui (x,t)) ) (6) 
P \OTi orj 



and 



Tij(r,t) = — (ui(x,t)uk(^,t)uj(x\t) - Ui(x,t)uk(^\t)uj(x',t)) (7) 

The tensor Ty comes form the inertia term in Navier-Stokes equation and, as it involves a product of third order in 
the velocity, reflects the fact that there is not a close set of equation for the correlations of successive orders but there 
is a hierarchy of equations instead. The problem of closing that hierarchy is known as the "moment closure problem" 
p2| . Let us call the Fourier transform of Rij(r,t). Then the energy density becomes 

^Rii{0,t) = J E(k,t) dk, 

where 

E(k,t) = ±J $ 2l (k,t) fc 2 dn(k) (8) 

is the energy density stored in eddies of size A: -1 . Defining as the Fourier transform of T^-, we obtain from (|5|) 
the equation of balance of the energy spectrum: 
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- -7rE(k, t) = T(k, t) + 2vk 2 E(k, t) 



(9) 



where 

T(k,t) = -±Jr u (k,t)k 2 dn(k) (io) 

The inertia term T(k,t) is the one that contains the mode-mode interaction, and its effect is to drain energy form 
the more energetic modes -typically the bigger ones- to the ones where there is major viscous dissipation -the smaller 
ones-. 



B. Flows in expanding universes 



Robertson - Walker (FRW) background 



dt 2 + a 2 {t) (dx 2 + dy 2 + dz 2 )), the generalization of the above arguments 



For a curved space-time, in particular that described by a Friedmann 
metric with zero spatial curvature (ds 2 

has been considered by many authors |46|-|50|] . We follow Tomita et aVs analysis |38|, in which they obtain the solution 
for the energy spectrum in the case of homogeneous, isotropic and incompressible turbulence. 

In a generic space time, we describe fluid flow from the energy density p, pressure p and four velocity U . The 
symmetries of the FRW solution suggest using instead the commoving three velocity u % — U l /U°; if U % <C U° the 



=0, it is incompressible (u = (try. 



x 3 )). 



Later on, we shall also use the physical 



flow is non relativistic, and if Vu 
three velocity v = a(t)u. 

The corresponding continuity and Navier-Stokes equations for a Robertson- Walker background are obtained by the 
condition of conservation of the energy-momentum tensor p(| . For a non relativistic incompressible fluid, with shear 
viscosity 77 = v (p + p) (but no bulk viscosity) , these reduce to: 



dp 
dt 



•3-(p + p) 
a 







(11) 



du 



(u-V) 



d In ((p + p)a 5 
dt 



u 



a 2 (p + p) a 2 



;A7 2 u 



(12) 



where we have assumed that p + p depends only on time. For the physical three velocity v, the corresponding 
Navier-Stokes equation reads: 



~dt 



1 enn((p + p)a 4 ) 

- (v-V) H '- 

a dt 



Vp 



1 



a (p + p) a? 



iA7 2 v 



(13) 



In obtaining (|ll|)-(|13|) we have neglected possible perturbations to the FRW metric. The corresponding equations 
considering fluctuations in the metric {g^ = g < t v + h IJ , v ) have been obtained by Weinberg J2(|. The continuity equation 
is not corrected by gravitational perturbations, while in the Navier-Stokes equation the metric fluctuations appear 
explicitly only within the shear viscosity term. It can be demonstrated that these terms involving metric fluctuations 
are negligible for scales that are inside the horizon (3^]. For scales bigger than the Hubble radius, since dissipation 
through viscosity is not effective anyway, we may still use the unperturbed Navier - Stokes equation. 

The operation of Fourier transforming in the case of a Robertson- Walker cosmology is done in terms of commoving 
wave-numbers. In doing so, the following equation for the energy spectrum is obtained: 




T(k,t)+2l — 



din ((jp + p)a 4 
dt 



E(k,t) 



(14) 



where the relationship between E(k,t) and 3>y(fc, t) as well as between T(k,t) and Tij(k ) t) is the same as that for a 
flat space time, if we define Rij and from correlations of physical quantities, as follows: 

R i3 (r,t) = a 2 (uj(x,£)uj(x + r,t)) , 



d 

Tij(r,t) = fl2 T^- ((«i(x,i)uk(x,t)iij(x + r,i)) - (ui(x,t)u k (x + r,t)uj(x + r,t))) 



(15) 
(16) 
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III. EQUIVALENT FLUID FOR INFLATON FLUCTUATIONS 



After establishing the basic necessary notions for the description of hydrodynamic flows, our goal is to associate an 
equivalent fluid description to inflaton fluctuations, and to derive the spectrum of primordial density fluctuations at 
reheating therefrom. We shall discuss in the following sections some non trivial instances of this method. 



A. The inflaton as a fluid 



To describe the inflaton field from the point of view of an equivalent fluid, we need to obtain the energy density, 
pressure and velocity of this fluid as functionals of the state of the field. To this end, our starting point will be that in 
the rest frame of the fluid (quantities in this frame being labelled by a curl), the field ought to be spatially constant 

V0 = O (17) 

To obtain the fluid four-velocity, we make a boost to the commoving frame. Then, the boost's characteristic velocity 
will be the one we are seeking for. By the condition (|l7]) we obtain: 

Ui = — (18) 



which is generalized to the covariant form 



«„ = - ; M (19) 

The energy density in the rest frame must be: 

Using the Lorentz transformations with the four-velocity ( fl9| ) we obtain the general form for the energy density: 

1 

P = —2 



d p cf> d p (f> + V{(/>) (20) 



Finally, we obtain the pressure imposing an equality between the energy-momentum tensor for a perfect fluid (see for 
example p6[ ) and that for a minimally coupled scalar field (see |l3|]). The resulting pressure is: 

p=-\d p <t>d?<t>-V{4>) (21) 

Since the possibility of deriving a Navier - Stokes equation for the equivalent fluid rests on the conservation of 
T^ v , in principle only the whole inflaton field can be thus represented. However, under the approximation that the 
homogeneous part of the inflaton essentially contributes an effective cosmological constant, the background energy 
momentum tensor Tq^ — Ag^ is independently conserved (even if A were not constant, conservation fails only on 
scales too large to be cosmologically relevant), and we can associate an equivalent fluid to the inhomogeneous quantum 
fluctuations 8(j> alone. For this fluid, we find the physical velocity (v l = au l — a^ui) 

vi=K hys (^j (22) 

where 4>o is the homogeneous background. We interpret this equation to mean that stochastic averages of the fluid 
velocity are to be identified with (symmetric) quantum expectation values of the operator in the right hand side 
[ pl| , [52"| |. As far as the equation of state is concerned, the free energy for a massive scalar field in the high temperature 
limit (T ^> m) |5j| gives us the relationship between the pressure and the energy density, which turns out to be 
that for radiation, p = This means that the energy density for this fluid redshifts proportional to a -4 . As this 
result has been obtained for a flat space-time, it is valid for scales smaller than the curvature radius. When scales 
are bigger than the Hubble radius, which in turn takes place when the high temperature limit is no longer valid, the 
field's equation of motion □</> + V ((f)) = turns out to be 



G 



^ +3 ^ +T/(0O)=O 

for the background field 4>o (where we have neglected spatial derivatives) and 

— -p h 3H— t 1" ^ (</>o)<><P = (23) 

for the fluctuation (50 (where we have neglected spatial derivatives as well as non linear terms). The time derivative 
of the equation for the background field gives us an equation of motion for O : 

^ + 3fl|+H^=0 (24) 

where we have used the constancy of H during inflation. Comparing (^3|) and (|24]) we see that, as (/>o and S<f> obey the 
same equation of motion, they must be related by: 5<f) — Sq f (r), which implies that the ratio ^ must be independent 

of time (cfr. [p4|). For our fluid description, this means that when scales are much bigger than the Hubble radius, 
the velocity Ui (equation (|l^)) must remain constant, which in turn means that the physical three velocity v l = au l 
must redshift proportional to a -1 . As these scales are frozen out because they are outside the horizon, they cannot 
interact among them or be dissipated by viscosity. Thus, the Navier - Stokes equation (|l3| ) reduces to 

gy ain(Q + p)<) 

dt dt K ' 

We obtain v cx a" 1 when (p + p) oc a -3 , corresponding to the equation of state of matter: p = 0. Thus, when the 
scales are well outside the horizon, our fluid behaves as pressureless dust, which agrees with the well known prediction 
based on Virial's theorem for the equation of state of a field undergoing oscillations filf , which occurs at the final 
period of inflation. We must point out that the hypothesis of incompressibility is no longer valid for an equation of 
state of this type. Nevertheless, for scales bigger than the Hubble radius, which cannot decay through non linear 
interaction or dissipation by viscosity, equation (|2^) is still valid, regardless of the ratio of typical velocities to the 
speed of sound. 

B. Transport coefficients 

The framework to obtain transport coefficients for our fluid is linear response theory. In the limit of slowly variations 
in space and time of the magnitudes involved in the equation of conservation for the energy-momentum tensor, the 
system's response while it is slightly displaced from equilibrium can be alternatively described by Navier-Stokes 
and continuity equations as well as by equilibrium expectation values of correlation functions. Matching these two 
descriptions, one obtains the Kubo formula for the shear viscosity [56): 



— lim 

6 w,k—>0 



- [dt [ d 3 v e l ( k - r -^<[^(r,i),^(0,0)]) 



eq 



(26) 



where 7Ty are the traceless spatial-spatial components of the energy-momentum tensor: 
For a minimally coupled scalar field the commutator involved in (09) turns out to be: 



[7ry(aO,7r«(0)] = [Mz)a^(a!),0V(OW(O)] - \ [^(x)^ <f>(x) , 0^(0)0^(0)] (27) 

These commutators can be evaluated from retarded and advanced Green functions for a massive scalar field coupled 
to other fields as well as to itself, which must satisfy: 
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where T is the thermal width which comes from the field's interactions. From these Green functions we obtain the 
Pauli- Jordan two point function G(x,x') = ([4>(x), 4>(x')]): 



G(X,X') = -—r / d 4 P 



2iY f e -i(po(t-t')-p-(x-x')) 



Po 



{2irf J " l(Po - iT) 2 - p 2 - m 2 ] [(p + iT) 2 - p 2 m 2 ] 

The mean value for the anticommutator of the field Gi(x, x') = ({<fi(x), 4>(x')}) is obtained through the Kubo-Martin- 
Schwinger's theorem p7j |: 

2iT f. e- i ( Po(t -*' ) - p - (x - x ' ) )p coth(^ /2) 



G 1 (x,x') = --—j / d 4 p 



(27rf J " [(Po - iT) 2 - p 2 - m 2 ] [(p„ + ^^) 2 - p 2 - m 2 ] 

^From these, using Wick's theorem [)58| and the c-number character of the mean value for the commutator, we obtain 
the commutator of the product of two fields as: 

([<t>(x)4>{x)^{x')4>{x')]) = ([ct>(x),<t>(x')}) mx),<t>{x')}) 

The first term involved in u2% turns then out to be 



dt d 3 r e l ^~ wt) ( [di^dj^x), 9*0(0)^X0)] ) 



_ 4r 2 , {p 2 (k-p) 2 + [p-(k-p)] 2 }( U ;-p o )p coth(/3p o /2) 



(27r) 5 J ' \(p — iT) 2 — p 2 — m 2 | 2 \(w — p — iT) 2 — p 2 — m 2 \ 2 

and a similar formula for the second term in (p7|). This gives for the shear viscosity in the high temperature limit 
(T>m,r): 

i] = const T V 2 (28) 

The thermal width T comes from the imaginary part of the self-energy. From dimensional analysis, it must be 
proportional to the temperature since the only relevant scale in the high temperature limit is the temperature itself 
(the constant of proportionality must be much less than unity for the consistency of the high temperature limit). 
Assuming that the only present interaction is the one coming from a u^ 4 term, V must include two a insertions, which 
means that T must be proportional to a 2 (we will estimate it as T ~ a 2 T). The shear dynamic viscosity becomes 

r, - (t 4 T 3 (29) 

A similar analysis allows us to evaluate the bulk viscosity, which turns out to be zero for a fluid with an equation of 
state of the type p = -|p in agreement with our previous assumptions. 



C. Conversion of hydrodynamic fluctuations into primordial density contrast 



Having described the quantum field as a fluid, we will analyze the resulting spectrum of density inhomogeneities. 
To do so, we assume that at some time t\ during the beginning of the inflationary phase, when all the scales relevant 
to cosmology were inside the horizon, the scalar field fluctuations were undergoing hydrodynamic fluctuations. Once 
the scales leave the Hubble radius, their energy cannot be dissipated by viscosity or by nonlinear coupling. Thus, 
equation (Oh means that they evolve according to: 



E{k 1 t> t out ) = E(k,t= tout) 



{{p + p)a*){t) 



(30) 



where the subscript "out" refers to the time when each scale leaves the horizon. 

The definition of E(k) (equation (||)) can be written in terms of the Fourier transform of the velocity; since the 
flow is statistically isotropic and homogeneous: 
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<«W(k')> = (f§) * 3 (k + k') 



(31) 



Combining (J30|) and (|3lj) we can obtain the r.m.s. value for the scalar field velocity at the time of reheating, which 
will be the r.m.s. value of the perturbation in the radiation's streaming velocity. This perturbation will in turn 
produce fluctuations in the energy density of radiation, which will evolve in the usual way. The theory of relativistic 
very large wavelength fluctuations predicts 6 ~ a 2 , where 5 = Sp/p is the density contrast, and thus 5 ~ H8, while 
the continuity equation yields i5 ~ w/( on a scale of physical size I. Consistency of these two pictures leads to the 
relationship between the velocity at the time of reheating and the fluctuation in the energy density as 



dp 
P 



^reh 



IK 



reh 



(32) 



where H re h is the Hubble parameter at reheating. Following this fluctuations up to the time they reenter the Hubble 
radius, assuming that their size is such that they are always unstable (they must be always bigger than the Jeans's 
length), they grow following the law (see for example Jl|,p6[): 



dp 
9 



reh 



^eq 
&reh 



(O-ent 
a eq 



tt2 ,2 $p 
n reh 1 ~T 



(33) 



reh 



where the subscript "eni" means the time each scale reenters the Hubble radius (the second equation holds even if 
the entering time occurs before matter-radiation equality). Combining equations ((|32|)-(|33"|)), we obtain the density 
contrast predicted by this theory at the time the modes reenter the Hubble radius: 



(SkSk> 



Hreh a rehE(k' \ ~ ^reh) 



5 3 (k + k') 



(34) 



This is the main result of this paper, as it relates the density contrast to a hydrodynamic variable. We shall see a 
nontrivial application of this formula in next section, but before, it is convenient that we pause to show explicitly how 
the familiar results relating to free field fluctuations are recovered in this language. Probably the most important 
feature of a theory where inflaton fluctuations are free is that each mode evolves independently of the other ones. 
Immediately after leaving the horizon they freeze, a situation that can be described phenomenologically assigning to 
the mode the effective equation of state of dust. This implies that 



E(k, t = t reh ) = E(k, t = t out ) 



a(t = t out ) 

a(treh) 



(35) 



ix x \ E(k.t — t out (k)) 3 , 

(<Weni = ^2 6 (k + k) 



(36) 



Let us compare this expression to the usual one in terms of quantum fluctuations. First we use Eq. (|3l|), neglecting 
any variation of H or of the velocities during reheating, to get 



(6 k S k 



< v i {k)v i (k'))\. = 



t=t out (k) 



(37) 



We now relate the physical velocity to field fluctuations according to Eq. (|22|); at t — t ou t(k), k p h ys — H, and this 
reduces to 



(SkSk' 



H 



\t=t out (k) 



(38) 



which is the conventional result |>4| . 

This shows the agreement between the fluid description and the conventional approach in this case, although of 
course it is only in the nonlinear case where we expect the hydrodynamic formalism to bring definite advantages. 
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IV. SELF SIMILAR FLOWS AND NONLINEAR FLUCTUATIONS 



In the previous sections we set up the general formalism whereby we can associate to the evolution of quantum 
fluctuations during inflation an equivalent fluid description, and derive the corresponding primordial density contrast 
from hydrodynamic variables. Of course, to put the formalism to actual use, we must be able to solve the Navier - 
Stokes equations, which is in itself almost as daunting as solving the fundamental quantum field theory. However, 
there is in the hydrodynamic case a century of lore to draw upon f30| , and some well tested approximations leading 
to relatively simple solutions. In this section, we shall demonstrate the equivalent fluid method by investigating the 
spectra resulting from one of these solutions, namely self similar flows. Towards the end of the section, we shall discuss 
the relevance of these solutions to actual cosmology. 



A. Self similar flows in flat and expanding universes 

As we have seen in the previous section (Eq.(|34|)), the key element in deriving the primordial density contrast is 
the energy spectrum E(k) (Eq.(||)), which is the solution of the balance equation (Eq.(^)). In it, the right hand side 
contains the viscous dissipation as well as the inertial force T(k, t). The overall effect of this term is to transfer energy 
from a given scale to smaller ones through mode - mode coupling; thus it is natural to model the action of the inertia 
term as a source of viscous dissipation, where the effective turbulent viscosity for a given mode depends on the motion 
of all smaller eddies pi] . By providing closure, that is, writing this effective viscosity in terms of the spectrum itself, 
a closed evolution equation for E(k) is obtained. Concretely, Heisenberg [53 proposed the ansatz 



I T(k', t) dk' = 2v{k, t) / E(k', t) k' 2 dk' 
o Jo 



(39) 



where 



u(k, t) = A J™ y — ^pf^ dk' (40) 

and A is a dimensionless constant. With this hypothesis (known as the Heisenberg hypothesis) as the solution to 
the closure problem, Chandrasekhar |3?J has obtained the energy spectrum for decaying turbulence, assuming that 
there is a stage in the decay where the bigger eddies have sufficient amount of energy to maintain an equilibrium 
distribution, thus requiring that the solution for the spectrum should be self-similar. With this consideration into 
account he obtained an energy spectrum: 

j5( M) = * M F (±fty (4i) 

V ' A 2 k 3 t 2 J t UovW V > 

where kg and to are initial conditions (namely, the wave number corresponding to the bigger eddy and its typical time 
of evolution). The function F obeys the equation 



, ,dF(x') 
F{x')-x' y ' 







f'OC 


dx' - < 


vk 2 t 4 


L 




(42) 



dx' 

which predicts a Kolmogorov type behavior for an inviscid fluid (R — > oo, R — A. ) in the ultraviolet limit: 

^ 

F(x) -> const x~ 5/3 (v = , x -> oo) 

While for nonzero viscosity: 

F(x) — > const x~ 7 (y ^ , x — > oo) 

In the infrared limit, F has the universal behavior F{x) = Ax (x <C 1), and thus we find a linear energy spectrum 
for k\ft <§; ko^/tQ. Chandrasekhar's self similar solutions are easily generalized to flows in expanding Universes. The 
dependence on time and wave-number for the self similar energy spectrum is J38{ 
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where the subscript i refers to "initial" and A and vt are respectively the Taylor's microscale and an average turbulent 
velocity, defined as: 

s 5 JTO^ -/**,*)* (44) 

Their time evolution must follow the law: 



The viscosity for our fluid, at least at high temperature, is given by Eq. (|29|). The equation which determines 
de function F(Xk) in ( ft3| ) turns out to be the same as in flat space time, Eq. (0), which means that assuming 
Heisenberg's hypothesis the spectrum is linear in k for scales much bigger than the Taylor's microscale ^|: 

E{k, t) = Avl ( ^ffif ) 2 A?fc for Afc « 1 (46) 



B. Nonlinear inflationary models 



We now want to place a self similar solution in the context of a inflationary scenario where, instead of regarding 
the inflaton fluctuations as free, we shall substitute them by an equivalent fluid, whose evolution we will assume to 
be self similar. We discuss at the end whether this last assumption is a reasonable one. 

As before, we will assume a duration of inflation close to the minimum value (iV m i n ~ 60, where N stands for the 
number of e-folds), which can be justified by the expected quadrupole anisotropy [p3| as well as by the ratio of the 
present to the critical density (see Appendix). By this assumption, a scale whose present size equals the horizon 
(~ 3000 Mpc) leaves the Hubble radius soon after the beginning of inflation. 

Unless in the free field case, here we cannot deal with each mode independently, but we must treat the whole flow 
subject to a phenomenological equation of state. Let us assume the self similar flow sets in at a time t\ when the 
temperature T » H (we discuss whether this is a suitable assumption below); and that the present horizon scale 
leaves the horizon at or around time t\. Then it is valid to use the high temperature limit for length scales close to 
the present horizon while they leave the Hubble radius during the inflationary phase. The fluid's equation of state in 
this limit is of the p = ^p type, which means that the product (p + p) a 4 remains constant throughout the universal 
expansion. The factor ({p + p) a i ) involved in ( |30| ) is then independent of the particular scale being considered 
within this group. Thus, by © and © we can obtain the energy spectrum for these scales while they are outside 
the horizon: 

E(k,t>t out (k)) = E(k,t = t out (k)) 
E (k ,t > t out (k )) E(k ,t = W(fco)) 

where (cfr. ©) 

E(k,t = t out )=Av^t 1 )X 2 (t 1 )k (48) 

A(ii) being the commoving Taylor's microscale at the time the self similar flow sets in. As at the initial time t\ the 
only relevant scale is the temperature, we expect the initial Taylor's microscale to be the inverse of the temperature 



2 We wish to point out an ambiguity concerning the meaning of Heisenberg's hypothesis in the case of curved spaces. For 
flat space time, the proportionality between the integral up to a certain wave number k of the inertia and the viscous forces 
is given by (B9h and (kid). In the case of a FRW space time, the autosimilar solution required by Tomita et at (^) needs a 
time dependent dimensionless constant A proportional to r/a 2 for the consistency of the solution. This product does remain 
constant only if the dynamic shear viscosity evolves in time proportional to a -2 . Thus unless this is the case, the solution we 
described looks like a natural curved space generalization of the Heisenberg - Chandrasekhar solution, but does not admit the 
same physical interpretation. 
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at that time, i.e. X p hys(ti) ~ T(h) • ^ ( |45| ) we can obtain the commoving Taylor's microscale at later times, such 
that most of the scales are still inside the horizon: 

A 2 (<) = A 2 (0 + 10 f \-^dt ~ A 2 (i x ) + 10-^M— — i-— (49) 

4 (p + p)a 2 (p + p)(ti)ifa 2 (ti) 



where we have used the fact that for a viscosity dependence upon time given by (|29|) the main contribution to the 
integral is given by its lower limit, provided there is a difference between the times t and t\ of more than an e-fold. 
Since 

v(h) 4 T(ti) 



a 



X 2 (t 1 )H(p + p)(h) H 

this means that the commoving Taylor's microscale freezes in its initial value, unless coupling is very strong. Let us 
take as reference scale ko in Eq. ( J47| ) the inverse Taylor's microscale (which is the last scale to leave the horizon in 
the high temperature regime). Then, from Eq. ( |30| ) 

(50) 



E(k ,t re h) ( a(Tph ys — H) 



E(k ,t = t out (k )) V a(t reh ) 



From Eqs. @ and (g), 



and so from Eq. Q) 



£(/,/,„,) ^ l " n ^ " : Y lrfU:^U:^~.. ( (51) 
a{t reh ) J \a{t reh )H J 



<W ent = ^^ 3 (k + k') (52) 



That is, a scale invariant Harrison - Zel 'dovich spectrum |3l|] with amplitude Vt- 

On the other hand, let us recall that the fluid velocity is related to the underlying field description by v ~ k p h ys 5(j)/(j). 
Now the typical wave number is k p h ys ~ T, and at high temperature also 6cj> ~ T. From the background evolution, 
we have <f> ~ V'{4>)/2>H ~ m 2 _ff/30. The number of e-foldings iV ~ H<P/4> ~ ( t )2 / m p: so ~ */Nm p , <f> ~ m p H, and 

- ^ (53) 

Since the ratio is necessarily less than one, the hypothesis of a non relativistic incompressible flow is seen to be 
consistent. The constraint of A; A -1 reduces to a minimum scale above which we obtain scale invariance. Jeans's 
length imposes a lower limit bigger than this (we are assuming that the scales are always unstable while they are 
outside the horizon, which is valid if they are bigger than the Jeans's length). 

Finally, combining the estimates for the initial Taylor's microscale and the turbulent velocity, we obtain a Reynolds 
number: 

R _ 4 \phys(ti)v t {t{)p(t{) _ Vt(h) 



3 r)(h) <J 4 

suggesting highly turbulent motion, specially for small couplings (the self similar solutions are not dependent on high 
Reynolds numbers anyway). The estimate of the viscosity must be taken with a grain of salt, nevertheless, since we 
have ignored non perturbative effects j5j| . 

Observations suggest that density fluctuations have a scale invariant spectrum with an amplitude Sp/p ~ 10 -5 over a 
range of scales going from maybe as low as lOOMpc up to the present horizon scale at l now = 3000Mpc ~ lO^GeV^ 1 
|jl5f |62) . Recall Eq. ( |57|) from the Appendix, relating the size of the Universe at the time t ex u when the actual 
horizon's scale left the inflationary horizon, to the Hubble parameter during inflation. At the time t ex u, the physical 
Taylor microscale is at most A = (a(t ex it) / a(tx)) T^ x (ti) = 10 _2 i? _1 , since a larger value would narrow too much the 
regime where the spectrum is scale invariant. So 

T(h) f a(t exit ) \ 1Q2 / H \ 1/2 = 1Q _ 26 f a{t f ) 
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We obtain agreement with observations (cfr.(p2) and (p3|)) provided 



«(*/) 



10 23 = e 53 



Since the total run of inflation is some sixty e- foldings (no more than 70 even in the extreme case H = m p ) this is 
not hard to obtain, provided self similarity sets in a few e-foldings after the beginning of inflation. 

We have thus set an explicit model where interacting inflaton fluctuations lead to a density contrast in agreement 
with observations. We could rest our case at this point, but before that, we would like to discuss briefly how likely it 
is that this solution was actually realized in the Early Universe. 



C. Self similar flows and our Universe 



As we have seen above, a self similar turbulent flow pattern in the equivalent fluid could explain the scale invariant 
spectrum of primordial density fluctuations observed at scales above lOOMpc, provided the self similar regime sets 
early enough. Our goal in this final section is to discuss whether this is a likely assumption regarding our own Universe. 

In general, it is known that turbulent flows tend to relax towards self similarity, but it is difficult to estimate how 
long does it take to get there. In a rough approximation, there are two issues involved, first, on what time scale the 
Heisenberg closure condition becomes valid, and then, how long does it take the closed equation for the spectrum to 
yield self similarity. 

If posed in these terms, we may observe the analogy with the problem of equilibration in the theory of dilute gases. 
In the latter, the fundamental description of the dynamics is the BBGKY hierarchy [Q. However, after a time of 
the order of a collision time the hierarchy can be closed, turning into the Boltzmann equation, which in turn leads to 
equilibrium in times of the order of the mean free time. 

In the turbulence problem, the time scale of a typical eddy is 



A _ m, p H 1 / H 



v t T 3 H \Tv t 



(54) 



We may take this as an estimate of the mean free time. The collision time can be estimated as the characteristic 
time for the smallest eddies in the flow (we visualize a collision as the exchange of a small eddy between two larger 
ones). According to Landau - Lifshitz fl64| , this is 

r rai „ = TiT 3 / 4 « r (55) 

It is simple to find parameters leading to r of the order of an e-folding, and therefore r m - m -C H \ for example, 
take H = 10 6 GeU, T{t x ) = 10 10 GeV, v t (h) = 1CT 5 , with T reh = 10 12 GeV. These solutions may not be available in 
models with such simple potentials as acj) , but they are feasible in multiparameter models such as those invoked in 
the potential reconstruction program p4| . While this does not constitute a proof, it makes it plausible that a self 
similar solution may be realized in the early stages of inflation. 

Once the motion becomes self-similar, it may last for a rather long time. From the equation of motion Eq. ( |l4| ) 
we expect the energy in the flow to be dissipated in a time scale of the order of r re i — v~ x }? phys . With A ~ T^ 1 and 
v ~ r\l (p + p) ~ <r 4 T _1 , we obtain 

1 . x 

Trel ~ ~ Rt (56) 

This time will be generally larger than an e-folding, and might be even larger than the whole duration of inflation 
if R is large enough. 

We conclude that self similarity will be easily achieved for scales around the Taylor microscale or smaller, and will 
propagate to larger scales after several e-foldings. At larger scales there could be deviations from scale invariance, 
associated with the transient behavior of the flow. We may conjecture that as the fluid cools down typical velocities 
will decrease, and thereby the primordial density contrast too. Thus the model will naturally yield higher power at 
larger scales, reproducing the behavior of multiple inflation models J65| . 

It bears mention that the same situation occurs in the usual treatments of the free field case. From Eq. (|3|), we 
may conclude that the spectrum of primordial density fluctuations has an amplitude 5k ~ (H 2 /<fi)\/l + njfe, where 
nk is the occupation number for that mode in the initial state of the inflaton field. Vacuum dominance requires 
rik <C 1 for k > ^Jm p H, but places no real restriction on larger scales which, in minimal models of inflation, are still 
observable. Thus we only obtain the conventional result (corresponding to n& = for all modes) for specially chosen 
initial conditions. 
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V. FINAL REMARKS 



In present conventional approaches, density inhomogeneities arise from primordial fluctuations in the inflaton field, 
ultimately of quantum origin. Fluctuations are treated as a free field, thus forcing upon us the assumption that higher 
derivatives of the inflaton potential are negligible. In this paper we sought a direct estimate of the primordial density 
contrast generated in a nonlinear inflationary model. Instead of assuming fluctuations to behave as a free field, we 
consider them to be coupled, so that they can be described phenomenologically as a fluid. We showed that there are 
flow patterns for this fluid that reproduce observations at very large scales. 

While our model takes advantage of many insights from earlier studies of the role of turbulence in structure formation 
(see |3l],|38| etc.), it explores a whole new aspect of the problem in the sense that it places turbulence at the origin 
of the primordial fluctuations, rather than being excited from them. It is therefore free from the criticisms that are 
usually raised against the turbulence theory of structure formation |2^,|3(J. To the contrary, our model successfully 
reproduces the results of the conventional approach on very large scales, that is, a scale invariant spectrum with a 
density contrast of about 10~ 5 . While it is not completely free of fine tuning, the fine tuning involved is the same as 
in the usual scenario. 

In our view, the main result of our work is not that a self similar solution should be our final description of 
fluctuations during inflation, but rather that it is possible to make sense of the physics of fluctuations even in rather 
general potentials. The self similar solutions we have explored in some detail should be seen as an ideal case which 
will more or less approximate actual flow patterns; indeed, the same could be said of the De Sitter invariant vacuum 
as a description of the actual state of the field in free theories. 

The connection of hydrodynamics to fluctuation generation has some interest of its own, as it provides an alternative 
to brute force quantum field theoretic calculations, and also yields physical insight on the macroscopic behavior of 
quantum fields in the Early universe. The equivalent fluid method may be used to advantage also in other regimes, 
such as the pre - inflationary Universe and the reheating era, where the strong back reaction phase has so far been 
untractable |^],[36| . Moreover, it opens up a wealth of new phenomena, such as intermittence |^3|,[54j and shocks [Q , 
which are not apparent in the customary treatments. We will continue our research in this field, which promises a 
most rewarding dialogue between cosmology, astrophysics, and nonlinear physics at large. 
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VII. APPENDIX 



The minimal amount of inflation necessary to solve the homogeneity problem is obtained by the condition that a 
scale of the size of the horizon at present (~ 3000Ajfpc) should have been inside the Hubble radius at the beginning 
of inflation. The Hubble radius during inflation is approximately constant. A scale whose physical size at present is 
X(to) was, at the end of reheating 

A(W) = A(to)%^ = A(t„)^ - A(t„) 2.35 ■ lfl-» f^T) 



Write T re h — Ty/m p H, where V < 1 and H is the Hubble parameter during inflation. Then towards the end of 
inflation we have 



A(t / )=A(* )2.35-10- 25 (^^y /2 
This scale left the horizon at a time t ex n with A (t ex u) = (a(t ex it) / a(t f))X (it) — H . So 



^h^Y] A(i )2.35 • W- 20 GeV ( 



H 



(t f ) J y VJ ■ \10 5 GeV 



1/2 
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In particular, for the present horizon scale we get 



a(t e 



*(*/) 



H 



10 5 GeF 



1/2 



(57) 



Therefore, defining iV m in — In [a(tf)/a(t ex i t )], which would make t ex i t — U for the scale of the horizon at present 

/ H \ 1/2 
V 10 5 Ge^ J 



= In 



= 64.4 



InW — 



(58) 



On the other hand, inflation could not have lasted much more than this because otherwise the present density should 
be so fine tuned to the critical one that would contradict observations as well as most speculations on dark matter's 
density. This can be seen from one of the Friedmann's equations: 



n(t) - 1 = 



k 



(a(t)H(t)Y 



where fl is the ratio of the density to the critical one and k is the spatial curvature of the FRW metric. This means 
that the ratio (assuming instantaneous reheating): 



Q(t ) 



1 



= (flfe)fffe)) 2 

n(*i)-l (a(t )H(t )f 



exp 2 



69.06 + In 




N 



(59) 



So, if we set £l(ti) to lie in the interval (0, 1), we conclude that if fl(to) is not fine tuned to 1, N should not have been 
much more than 60. 

There is another argument supporting that inflation should not have spanned much more than its minimum duration 
based on the comparison between the expected quadrupole anisotropy and the detected one (see [p3| ). 
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